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I.  INTRODUCTION 


Wide-angle  biconical  antennas  have  been  used  for  about  two  decades  to  produce 
both  vertically  and  horizontally  polarized  electromagnetic  fields  for  use  in  electromagnetic 
pulse  (EMP)  simulation.  Approximate  theoretical  models  have  been  used  in  conjunction 
with  field  mapping  to  provide  a  practical  basis  for  testing  and  evaluating  military  systems 
against  the  EMP  threat. 

In  recent  years,  however,  there  has  been  a  greater  interest  in  the  development  of 
biconical  antenna  EMP  simulators  that  can  deliver  pulses  with  very  fast  risetimes.  A  major 
element  of  the  waveform  design  for  this  type  of  simulator  is  the  resistive  loading  of  the 
antenna,  which  predominantly  affects  the  intermediate-  to  late-time  character  of  the  radiated 
field.  Over  the  years,  numerous  models  have  been  developed  to  approximate  the  effects  of 
resistive  loading  on  field  behavior.  However,  to  my  knowledge,  this  investigation  is  the 
first  attempt  to  develop  a  rigorous  theory  leading  to  the  prediction  of  the  far  fields  from  a 
resistively  loaded  wide-angle  biconical  antenna. 

Schelkunoff  [1]  appears  to  have  been  the  first  to  develop  an  exact  theory  of  the 
biconical  antenna  in  the  absence  of  resistive  loading.  His  method,  although  rigorously 
correct,  presented  formidable  computational  difficulties  at  the  time  of  its  inception  (1941). 
The  principal  difficulties  appeared  to  be  the  computation  of  the  0-dependent  eigenfunctions 
that  are  characteristic  of  the  biconical  antenna  with  perfectly  conducting  interior  surfaces. 
As  shown  in  the  appendix,  this  calculation  requires  the  determination  of  the  roots  of 
polynomial  equations  involving  hundreds  of  terms.  While  this  computation  is  feasible 
today,  it  was  virtually  insurmountable  in  Schelkunoff  s  time. 

Because  of  these  numerical  problems,  emphasis  was  placed  on  obtaining 
approximate  solutions  in  selected  regimes.  For  example,  Tai  [2]  developed  a  method  for 
small-angle  bicones,  which  is  also  briefly  summarized  by  Kraus  [3].  Techniques  relevant 
to  wide-angle  bicones  were  developed  by  Smith  [4]  and  Tai  [5].  In  all  these  cases, 
however,  the  issue  of  resistive  loading  was  not  addressed;  perfectly  interior  conducting 
boundaries  were  assumed. 

The  purpose  of  this  investigation  is  to  theoretically  evaluate  the  possibility  of 
rigorously  solving  the  biconical  antenna  problem  with  resistive  loading.  It  appears  that,  at 
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least  in  a  formal  sense,  I  have  been  successful  inasmuch  as  the  methodology  leads  to 
closure  on  a  solution.  My  method  is  an  extension  of  the  Schelkunoff  formalism  in  the 
absence  of  resistive  loading.  Since  this  case  forms  the  basic  building  block  for  the  method, 
I  provide  a  comprehensive  discussion  of  Schelkunoff’s  methodology  [6]  using  Kong’s 
formalism  [7].  This  review  is  rendered  in  section  2. 

In  section  3  I  extend  the  theory  for  the  case  of  a  resistively  loaded  bicone.  The 
resistively  loaded  case  is  found  to  require  additional  terms  in  the  antenna  region,  compared 
to  the  lossless  case,  so  that  the  interior  boundary  condition  is  satisfied.  The  resulting 
equations  do  not  lend  themselves  to  analytical  solutions  in  closed  form.  Closure  is 
achieved  through  the  introduction  of  an  orthogonal  basis  of  radial  functions,  which 
ultimately  reduces  the  problem  to  one  that  can  be  solved  in  matrix  form.  The  computational 
implementation  of  this  theory  remains  to  be  developed. 

A  key  element  in  the  computational  algorithm  is  the  evaluation  of  the  0-dependent 

eigenvalues,  m/,  and  eigenfunctions,  (cos  do)  of  the  biconical  antenna  in  the  absence  of 

i 

resistive  loading.  These  entities  are  essential  in  determining  the  parameters  of  the  matrix 
equation  leading  to  the  determination  of  the  fields.  The  computational  considerations  for 

(cos  6lo)  are  presented  in  the  appendix. 
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2.  REVIEW  OF  SCHELKUNOFF’S  METHOD  USING 
KONG’S  FORMALISM 


This  section  establishes  the  basic  formalism  that  will  be  used  in  the  general 
treatment  of  the  resistively  loaded  biconical  antenna  of  section  3.  This  material  is  extracted 
from  Kong’s  text  [7].  No  attempt  is  made  to  reproduce  the  step-by-step  derivation  of  the 
results.  The  reader  is  referred  to  Kong’s  text  for  more  details.  Wherever  possible  I  have 
used  Kong’s  notation,  the  notable  exception  being  the  use  of  for  his 

This  summary  is  rendered  in  a  manner  which  highlights  those  changes  that  are 
necessary  to  extend  the  theory  from  the  lossless  case  to  the  resistively  loaded  bicone. 
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Figure  1  shows  a  model  of  the  biconical  antenna,  and  the  corresponding  coordinate 
system.  For  the  cases  of  interest  only  radial  currents  are  present,  with  the  accompanying 

conditions 


Figure  1.  Geometric  considerations  for  biconical  antenna. 


Under  the  foregoing  conditions.  Maxwell’s  equations  reduce  to 


-ll] ■ 

(3) 

(4) 

(5) 

where  we  have  replaced  8/Bt  by  jco. 

In  the  air  region,  defined  by  the  regime 

r>L  =  radius  of  bicone  , 

the  solution  of  equations  (3)  to  (5)  using  separation  of  variables  yields 

N=1 

iZ,  (2)  • 

^  ^N(N  +  DVa,  (Wa,(cos  ff)  ,  (S) 

N=l 

N=l 

(2) 

where  Psicos  0)  is  the  Legendre  polynomial  of  order  A^,  (kr)  is  the  Hankel  function  of 
the  second  kind  of  degree  N,  the  bi^’s  are  constants  to  be  determined  from  the  solution  of 

the  problem,  Zq  is  the  impedance  of  free  space  =  and 

Pj^  (cos  0)  =  dP/y  (cos  0)/d0  .  (10) 

It  is  important  to  note  that  the  functional  form  of  equations  (7)  to  (9)  will  be  the 
same  with  or  without  resistive  loading.  The  mathematical  structure  of  the  external  fields 
stems  from  the  requirement  that  only  outgoing  waves  be  present  outside  the  antenna,  as 
well  as  the  condition  of  symmetry. 


(6) 

(7) 
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H(7c-e)  =  H(e)  , 


E^{n-d)  =  E^{e) 


When  equation  (11)  applies,  only  odd  values  of  are  allowed  in  the  summations  of 

(2) 

equations  (7)  to  (9).  The  first  Hankel  function,  ,  is  given  by 


1-^  ^ 


kr\  kr  • 


while  the  first  Legendre  polynomial  and  its  derivative  are  given  by 

Pi  =  cos  d  ,  (13) 

Pj  =  -sin  e  .  (14) 

Using  equations  (12)  to  (14),  the  leading  terms  of  Eg,  and  Er  in  the  far  field 
{kr  »  1)  become 

=  (15) 


•  (17) 

{kr) 

Examination  of  equations  (15)  to  (17)  shows  that  for  practical  purposes  a  determination  of 
b\  will  be  sufficient  to  determine  the  fields  in  the  test  volume. 

After  some  mathematical  manipulation  it  can  be  shown  that  the  interior  fields,  valid 
in  the  region  ^  <  0  <  ;r-  can  be  written  in  the  form 

=  2^  *  U  S  e)  ,  (18) 
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(19) 


A. 

iTckr 


^  Uiu+l)  aJJJcr)TJf:os  B)  , 


where 


Vo  (^) 

Z^lnr  sin  6 


+ 


A 

2;0i:r 


7’u  =  ^  ^  ~  /*u(-cos  0)]  , 


(20) 


(21) 


=  dJJdd  .  (22) 

V  (L) 

/o(/-)  =  [j  sin  k(L  -  r)  +  cos  k(L  -  r)]  .  (23) 

(r)  =  V(j  (L)  [cos  k{L  -  r)  +jY^,  sine  kit  -  r)j  ,  (24) 

Zq  In  cot 

Zc  = - - — —  =  characteristic  imjjedance  of  bicone  ,  (25) 


juilif)  is  the  spherical  Bessel  function,  and  the  OuS  are  constants,  which,  like  the  bs's  of 
the  exterior  region,  are  to  be  determined  from  the  solution  to  the  problem.  The  terminating 
admittance  Yt  is  likewise  determined  from  the  solution. 

In  contrast  to  equations  (7)  to  (9),  in  which  the  summation  index  is  defined,  the 
values  of  u  in  equations  ( 1 8)  to  (20)  are  determined  from  the  boundary  condition  in  the 
antenna  region.  In  the  absence  of  resistive  loading  we  require  that 

Er{Bo)  =  Er{n-  do)  =  0  .  (26) 

Based  on  equation  (19),  the  specific  values  of  u  which  will  satisfy  equation  (26)  are 
solutions  of  the  equation 

r„(cos  Oq  )  =  T„(cos  {n-  ))  =  0  , 

which  reduce  to 

P„(cos  0Q )  -  P„(-cos  ^ )  =  0  . 


(a) 


(b) 


(27) 
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The  Pu’s  of  equation  (27)  are  not  Legendre  polynomials,  but  are  Legendre 
functions  since  u  is  not  generally  an  integer.  Appropriate  analytical  expressions  for  Pu  that 
are  necessary  to  solve  equation  (27),  and  hence  determine  the  values  of  u,  can  be  found  in 
Schelkunoff  [6],  Abramowitz  and  Stegor  [8],  and  Erdelyi  et  al  [9].  For  the  time  being  we 
simply  assume  that  equation  (27)  has  been  solved  and  the  values  of  u  determined.  Thus, 
the  summation  index  of  equations  ( 1 8)  to  (20)  is  defined. 

The  complete  solution  to  the  problem  is  now  found  by  matching  the  boundary 
conditions  at  r  =  L.  Using  equation  (23)  we  have 

I^(L)  =  Y,V^(L)  ,  (28) 

which  when  inserted  in  equation  (18)  gives 

YV^(L)  ,  ^  , 

H(r  =L)  =  — ^  >  a,j  ,  {kL)  T  .(cos  6)  .  (29) 

^  2;cLsin  0  2;r  ^  ^  ^ 

U 

Multiplying  both  sides  of  equation  (29)  by  sin  6  7^,  (cos  0)  and  integrating  from 
dQto  jT-  Oq  gives 

n  -6l(j 

"  TrjfkL)  ’  (30) 

where  we  have  used  the  normalization 

sin  0T„(cos  0)T  .(cos  0)  =  I  ^  (31) 

“  lN^ifu  =  u'  ^  ^ 

Equation  (31)  is  a  general  property  of  the  Legendre  functions  which  satisfies  the  boundary 
condition  of  equation  (27),  with  Nu  being  the  associated  normalization  constant. 

Equation  (30)  provides  a  connection  between  the  Ou’s  and  the  b^’s  through 
substitution  of  equation  (7).  We  have 

oo 

'  (32) 

N=\ 


i 


where 


a 


uN 


\  ikL) 
Nj,(kL) 


sin  6  T„(cos  9)  Pf^icos  6)  dd  . 


(33) 
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We  should  also  recall  in  passing  that  the  functional  form  of  and  Eq,  as  given  by 
equations  (18)  and  (20),  respectively,  ensures  that  for  both  exterior  and  interior  regions  the 
following  conditions  will  be  satisfied: 

H^{6)  =  H^n-d)  ,  (a) 

Ee(e)  =  E0(n-e)  .  (b)  (34) 

A  ^ 

If  we  now  define  Ou  and  to  be  the  components  of  column  vectors  a  and  b , 
respectively,  and  otyv  to  be  the  elements  of  a  matrix  A,  we  can  write  equation  (32)  in  the 
compact  matrix  form 

a=Ab  .  (35) 

The  second  relationship  between  the  a^’s  and  b^'s  is  obtained  by  requiring 
continuity  of  Eq  at  the  boundary.  From  equation  (9)  we  have 

LEJ^r  =  X  b^G^{kL)P\{co%  &)  ,  (36) 

^  N=\ 


where 


(kL)  = 


(2) 

rh^  (kr) 


r-L 


(37) 


Multiplying  both  sides  of  equation  (36)  by  sin  6P^{cos  0)  and  integrating  from  0  to  ;r 
gives 


b 


N 


2NiN+  1) 
2N+  1 


Ink 

jZ^Gf^ikL) 


f. 


sin  OLEgir  =  L)  P^^icos  6)  dO 


where  we  have  used  the  relationship 


(38) 


r.2N-H 

J  sin  edPP^icos  6)  P^.icos  0)  =  i  +  1) 


ifiV  =  N' 


Using  equation  (20)  applied  for  r  -  L  in  equation  (38)  yields 

Vod^) 


(39) 


(40) 
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(41) 


where  the  relevant  constants  are  defined  as  follows: 


g  _2A^(A^+  1)  FJ,kL) 
Pnu  2N+1  Gj^ikL) 


sin  d  Pj^icos  0)  r„(cos  0)  d0 


^  _  ,AN{N+l) 
N-  J  2N+1 


k 

G^m 


Pj^icos  0o)  . 


(42) 


FuikL) 


(43) 


The  derivation  of  equation  (40)  was  based  on  the  requirement  that  E0{r  =  L)  =  O  when  0 
lies  outside  the  limits  0Q<0<n-0o. 

.  A 

By  identifying  as  the  component  of  the  column  vector  K  and  pNu  as  the 
element  of  the  matrix  B,  we  can  recast  equation  (40)  in  the  form 

b  =  Ba+(yQ(L)lZ^)K  .  (44) 


The  termination  impedance,  Yi ,  can  now  be  determined  from  equations  (35),  (44), 
(29),  and  (7).  This  is  accomplished  by  first  integrating  equation  (29)  from  Ooto  n-  Oq- 
We  have 


‘  4Vo(^^) 


H  (r  =  L)  de  . 


(45) 


We  now  substitute  equation  (7)  into  equation  (45)  to  obtain 

=  VvW^(cos0,)  .  (46) 

o'  ^=1 

The  foregoing  expression  can  also  be  cast  in  matrix  form  by  introducing  the 

^  (2) 

transpose  of  the  column  vector  ^  whose  components  are  ikL)Pj^{cos  0q)  .  Thus,  we 
write  equation  (46)  in  the  form 

y  _  _L 

‘  Z^q(L)  ^  ^  ’  (47) 


where  ^  is  the  transpose  of  i  and  is  a  row  vector. 
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In  an  actual  computation  we  will  truncate  the  series  for  u  and  the  series  for  iV  after  a 
finite  number  of  terms.  However,  the  maximum  number  of  terms  used  in  the  respective 
series  may  differ  since  the  number  of  terms  required  to  describe  the  fields  in  the  interior  and 
exterior  regions  may  not  be  the  same.  This  does  not  place  any  restriction  on  the  theory  or 
the  method  of  solution.  It  is  convenient  in  this  report  to  conceptually  regard  the  u  and  N 
series  as  being  truncated  after  the  same  number  of  terms,  say  n.  Let  us  assume  that  this  is 


done  so  that  A  and  B  are  now  iixn  matrices  and  a,  b,  K,  and  |  are  n-dimensional  vectors. 
We  now  substitute  equation  (35)  into  equation  (44)  to  obtain 


Ub  = 

Z. 


(48) 


where  U  is  the  matrix 


U  =  T-BA 


and  I  is  the  identity  matrix.  The  solution  of  equation  (48)  is 


-1 


(49) 


(50) 


where  U  Ms  the  inverse  of  U . 

Substituting  equation  (50)  into  equation  (47)  gives  the  following  expression  for  the 
admittance: 


•7-2 

Zf.1t 


(51) 


Examination  of  the  terms  of  equation  (51)  shows  that  Yt  is  a  function  of  the  bicone  angle 
Bo,  radial  dimension  L,  and  wavenumber  (frequency)  k.  For  specified  values  of  L  and  6b, 
the  admittance  can  be  expressed  as  a  function  of  frequency  co. 

Once  Ytico)  is  determined  from  equation  (51),  the  remaining  calculation  proceeds  as 
follows:  We  initially  determine  Vo(L)  from  a  knowledge  of  the  source  voltage  Vo(0). 
Using  equation  (24)  we  have 

Vq  (L)  =  (cos  kL  +  jYf  Zf.  sin  klf^Q  (0)  .  (52) 

Substituting  equation  (52)  into  equation  (50)  gives  6 ,  and  the  external  fields  are  then 
determined  from  equation  (7). 
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3.  SOLUTION  WITH  RESISTIVE  LOADING 


The  deduction  of  the  solution  with  resistive  loading  is  not  an  obvious  extension  of 
the  lossless  case.  A  fundamental  aspect  of  the  lossless  case  was  that  the  interior  fields  were 
made  up  of  modes  with  index  u  determined  from  the  solution  of  the  interior  boundary 
condition 

£,(0o)  =  £r(^-6b)=O  .  (53) 

This  led  to  equation  (27)  and  the  determination  of  the  values  of  u. 

When  resistive  loading  is  considered,  the  new  boundary  condition  at  the  walls 
becomes 

E,(do)  =  E,i7t  -  0o)  =  2 Jr)  /  Jr)  (54) 

where  /  Jr)  is  the  wall  current  and  2  Jr)  is  the  surface  impedance.  The  wall  current  is 
given  by 

/Jr)  =  27cr  sin  (^o)  •  (55) 

The  existence  of  resistive  losses  requires  modification  of  the  functional  form  of  the 
fields  in  the  interior  region.  In  the  absence  of  resistive  loading,  the  interior  fields  are  given 
by  equations  (18)  to  (20)  with  the  assumed  noninteger  values  of  u  being  determined  from 
the  solution  of  equation  (27). 

When  the  boundary  condition  of  equation  (54)  replaces  that  of  equation  (53) 
because  of  resistive  loading,  the  structure  of  the  fields  given  by  equations  (18)  to  (20)  does 
not  appear  to  be  sufficient  in  itself  to  solve  the  problem.  I  have  not  been  able  to  find  a 
means  to  uniquely  determine  the  values  of  u  that  satisfy  equation  (54)  and  the  other 
boundary  conditions  using  the  field  representation  of  equations  (18)  to  (20). 

It  appears,  however,  that  a  unique  solution  to  the  problem  can  be  obtained  if  one 
adds  to  the  field  components  of  equations  (18)  to  (20)  additional  terms  involving  integer 
values  of  the  summation  index,  which  are  also  solutions  of  the  basic  equations  in  the 
antenna  region.  These  terms  are  necessary  to  satisfy  the  boundary  condition  of  equation 
(54)  and  do  not  exist  in  the  lossless  case. 

The  deduction  of  the  0  dependence  of  equations  (18)  to  (20)  was  based  on  the  fact 
that  the  two  linearly  independent  solutions  for  //^  in  the  6  dimension  are  derivatives  of 
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Legendre  functions  of  the  first  and  second  kind,  Py(cos  6)  and  QuCcos  ff),  respectively 
[6-9].  In  the  condition  where 


t>  =  u  integer  , 

(56) 

the  function  Qu  satisfies  the  equation 

Qu(cos  6)  =  PJ-cos  8)  . 

(57) 

Using  the  result 

cos(;r  -&)=  -cos  6  , 

(58) 

combined  with  equation  (57)  and  the  symmetry  requirement 

(59) 

led  to  the  choice  of 

Tu  =  j(Pu-Qu>=j  [Pu  (cos  d)-P^  (-cos  e )) 

(60) 

as  the  only  candidate  0-dependent  solution  (compare  equation  (21)). 

When  we  allow  d  to  be  an  integer,  the  0-dependent  part  of  the  solution  changes  as 

follows.  For  the  case  where 

v  =  N  =  integer  , 

we  have 

P^(cos  ^  (^  +  9)!  ^ 

1^0  iN-q)\{q\r  2 

(61) 

N  p  p 

g^(cos  8)  =  P^(cos  0)  ln(cot^)  -  ^  . 

m=l 

(62) 

Using  the  foregoing  equations  we  easily  see  that 

P^(-cos  0)  =  (-l)^P^(cos  0)  , 

(63) 

Qj^(-cos  0)  =  (-1)^'^^  G^(cos  0)  . 

(64) 

Employing  the  symmetry  condition  of  equation  (59)  requires  that  only  odd  powers 
of  A/  be  used  for  the  Pj^  terms  and  even  powers  of  N  for  the  Qj\/  terms. 
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When  integer  terms  are  added  to  equations  (18)  to  (20),  the  expressions  for  the 
fields  can  be  written  in  the  form 

u  nt 

Er  =  ^ u (w  +  1) aj^ ikr) T^(cos  6)  +  ^ m(m  +  1)  aJ„ikr)rn(cos  $)  ,(66) 

u  m 


r^(cos  0)  =  P„icos  0)  if  m  =  odd  ,  (a) 

r^(cos0)  =  Q^(cos  0)  if  m  =  even  ,  (b)  (70) 

and  the  d^'s  are  additional  constants  to  be  determined  from  the  problem  in  the  resistively 
loaded  case.  The  index  m  is  an  integer. 

We  assume  that  at  the  boundary  6=  Oq.  In  addition,  we  continue  to  retain  the 
equation 

r„(cos^)  =  0,  (71) 

so  that  the  noninteger  values  of  u  remain  unchanged,  and  hence  are  known. 

The  relationship  between  the  d^'s  and  d^'s  is  determined  from  the  boundary 
condition  of  equation  (54).  Using  equations  (23)  and  (55)  we  have 

y  V  (L)Z“‘  sin  kiL-r)  +  V^{L)Y,  cos  k{L  -  r) 
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(72) 


+  Z^r sin  Oq  aj^ikr) f^(cos  9^ )  +  Z^r sin  6{j  Y  t,(cos  B^)  = 

u  J  m 

^  m(m  +  1)  aJJJcr)Xn  (cos  B^ )  . 

m 

There  does  not  appear  to  be  a  simple  way  to  relate  the  a«’s  to  the  a«’s,  but  it  can  be 
accomplished  in  matrix  form  by  assuming  that  in  the  region 

0<r<L  ,  (73) 

we  can  expand  all  the  r-dependent  functions  in  an  orthonormal  basis  using  a  complete  (and 
as  yet  unspecified)  set  of  functions  0n(r)  which  satisfy  the  standard  conditions 

«r)«r)<ir=5^  ,  (74) 

Jq 

where  n  and  m  are  integers  and  dmn  is  the  Kronecker  delta. 

The  aforementioned  procedure  is  accomplished  by  (1)  multiplying  both  sides  of 
equation  (72)  by  r,  (2)  setting 

y{r)=ZJZo  , 

and  (3)  using  the  following  relationships: 

oo 

rysink(L-r)  =  '^  d^r)  ,  (75) 

rt=0 

oo 

rrcos  k(L  -  r)  =  ^  fMr)  ,  (76) 

n=0 

oo 

y'^y„(*r)  =  ^  gnu<Pnir)  ,  (77) 

n=0 

oo 

(^'•)  =  X  ’  (78) 

rt=0 
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(79) 


■go 


where  the  dnS,f„'s,  gnu’s,  gnm’s,  and  ^nm’s  are  constants  determined  from  the 
orthogonality  condition  of  equation  (74).  For  example, 

gmn=  {  •^Jmikr)(l>nir)dr  .  (80) 

•'0 

Since  the  0n’s  form  a  complete  set,  the  substitution  of  equations  (75)  to  (79)  into 
equation  (72)  yields 

jVQ(L)Z~^  d„  +  Vq(L)YJ^  +  sin  ©qX  grufu^^os  ^)a„ 


+  sin  00  X  ^  X  •  (81) 

ffl  Ifl 

By  defining  dn,/n,  flu,  and  ^to  be  components  in  the  same  dimensional  vector 
space,  we  can  cast  equation  (81)  in  the  matrix  fOTm 


A  =  A  =  A 

J*  -  ^  W  W 


d+VQ(l.)YJ-^Ga  +  G  a  a  , 


where 


^nu  sin  Oo^nuT’ufcos  0o)  . 


Gnm  =  sin  rjcOS  ©q  )  . 


^ nm  2,jck  i)^wn  > 


The  solution  of  equation  (82)  which  renders  a  in  terms  of  a  is  given  by 


d=ya^jV^{L)Z,X  +  V^{L)Y^T]  , 


where 
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v  =  (h-g)  , 

(85) 

X  =  {H-§fd  , 

(86) 

(87) 

Using  equation  (84),  the  determination  of  the  admittance  Yt  and  the  constants 

a,  a ,  and  b  follows  in  a  manner  analogous  to  the  lossless  case  considered  in  the  previous 
section  (compare  the  analysis  beginning  with  eq  (28)).  For  brevity,  some  of  the  obvious 
intermediate  steps  will  be  omitted.  From  equation  (65)  we  have 


H^(r  =  L)  = 


y,Vo(i^) 

ItcL  sin  6  ^  In 


hjt  X  (cos  ^  ^ X  in (^=05  &)  (88) 


If  we  now  multiply  equation  (88)  by  sin  6  (cos  6)  and  integrate  from  do  to 
n-  do,  and  then  use  equation  (7)  for  //^(r  =  L),  we  obtain  the  following  result: 


A  A  A 

a  =  Ab  -Q  a 


(89) 


where  A  is  the  previously  defined  matrix,  and  fi  is  a  matrix  whose  components  are 


Q  = 

^  Nj^m 


M  V 

j  sin  0T^(cos  d)  (cos  d)  dd 


(90) 


From  equation  (38)  we  have  the  relationship  which  relates  b^  to  the  interior  fields. 
However,  we  must  now  use  equation  (67)  for  Eq.  Inserting  equation  (67)  into  (38)  gives 

~  X  ^  ^  >  (91) 

u  m 


where 


B  _2N(2A^+1)  F„{kL) 
2N  +  1  G^ikL) 


ji-80  ^  ^ 

I  sin  dPf^  (cos  6)r;„(cos  d)  dd  . 
fib 


(92) 
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It  is  observed  that  equation  (91)  is  similar  in  mathematical  structure  to  equation 
(40);  the  matrix  equivalent  becomes 

b=Ba  +  K^VoiL) Z~\  B  J  ,  (93) 


which  is  likewise  analogous  to  equation  (44).  B  is  the  matrix  whose  elements  are  given  by 
equation  (92). 


The  final  required  equation  for  determining  the  unknowns  in  the  system  is  found  by 
integrating  equation  (88)  between  6b  and  jc—Oq.  There  results 


Z,Vo(L)  ^  ^  Z,Vo(L)  n 


(94) 


where  V'  is  a  column  vector  whose  components  are  /„,(kL)r^(cos  6^). 

A  ^  - 

In  summary,  there  are  four  unknowns  in  the  problem:  Y^,  a  ,b ,  and  a .  These  are 
determined  from  equations  (84),  (89),  (93),  and  (94).  The  simultaneous  solution  of  these 
matrix  equations  involves  many  intermediate  steps  which  for  brevity  are  not  presented. 
The  final  result  is 

-afx^  +av%  +jaz;'fx, 

= -  ,T-  .  ;r - ^  •  <«) 

l  +  a^  AfjXj-av^Xj 

where 

Z^L 

a  =-|-  .  (a) 

X^  =  i!  -BAf^kz,  ,  (b) 

X2=AXj  ,  (c) 

M^  =  (r-BA)~^  (B-BQ)  ,  (d) 

M^=AM^-a  ,  (e) 

X^  =  if-VM^)~^VX^  ,  (0 
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A  A  ^  A  —1  ^  A 

X^=X^+M^X^+jZ,M^X^  . 


(g) 

(h) 

(i)  (96) 


.A 

Once  Fj  is  computed,  the  other  parameters,  a,  6,  and  a,  are  detemuned  from 
equations  (84),  (89),  (90),  and  (94). 


4.  CONCLUSION 


In  this  study  we  have  demonstrated  that  a  methodology  for  calculating  the  fields  of 
a  biconical  antenna  with  resistive  loading  is  theoretically  and  computationally  feasible.  The 
technique  draws  on  the  modal  analysis  concept  originally  developed  by  Schelkunoff  for  the 
case  without  resistive  loading.  However,  the  extension  to  the  resistively  loaded  case  is 
considerably  more  complex  from  both  a  conceptual  and  computational  viewpoint. 

The  resistively  loaded  case  is  found  to  require  the  existence  of  additional  terms  in 
the  antenna  region  (compared  to  the  lossless  case)  so  that  the  interior  boundary  condition  is 
satisfied.  The  solution  of  the  interior  boundary  value  equation  appears  to  require  the 
introduction  of  an  orthogonal  basis  of  functions,  0/i(r),  defined  in  the  range  0  <  r  <  L, 
where  r  is  the  radial  coordinate  and  L  is  the  bicone  radius.  We  have  not  as  yet  selected  the 
0n(r);  however,  this  does  not  limit  the  theoretical  analysis.  Using  the  basis  functions, 
0n(O.  a  matrix  formulation  is  developed.  In  order  to  achieve  a  practical  solution  to  the 
problem,  it  will  be  necessary  to  assume  a  finite  number  of  terms.  This  number  has  not 
been  determined,  but  will  surely  depend  on  the  bicone  angle  6o  and  radius  L,  and  the 
magnitude  and  radial  distribution  of  the  resistive  loading. 

In  summary,  the  implementation  of  the  technique  developed  in  this  investigation 
requires  the  selection  of  0n(^)  and  the  evaluation  of  certain  mathematical  functions  and 
integrals  which  depend  on  0/,(r).  When  this  is  accomplished  it  should  be  a  relatively  easy 
matter  to  predict  the  near  and  far  fields  from  a  resistively  loaded  biconical  antenna. 
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APPENDIX: 


DETERMINATION  OF  THE  ROOTS  m,  AND 
EIGENFUNCTIONS  r„,. 
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APPENDIX:  DETERMINATION  OF  THE  ROOTS  Ui  AND 

EIGENFUNCTIONS  Tu, 


In  this  appendix  I  address  the  determination  of  the  roots  of  the  equation 

(cos  Oq)  =4-  (cos  do)-Pu  (-cos  ^o)]  =  (^  • 

1  Li  i 


(A-l) 


where  is  the  Legendre  function  of  order  Wj,  and  6b  is  the  bicone  angle  of  figure  1  in  the 
body  of  the  report.  I  also  discuss  the  determination  of  the  6-dependence  of  ':he  function 


T  (cos  6)  =  4  [Pu  (cos  9)  -  (-cos  6)]  . 

<  Li  • 


(A-2) 


Letting 


we  write  equation  (A- 1 )  in  the  form 


zo  =  cos  flo  , 


(A-3) 


(A-4) 


Since  u,  is  not  generallly  an  integer,  the  complete  series  expansion  for  Pu  must  be  used. 
This  series  can  be  deduced  from  the  hypergeometric  function  [1]  through  the  relationship 


Pu(^o)  =  w  +  1;  1 ; 


(A-5) 


which  is  valid  in  the  range  I  1  -  zqI  <  2.  This  latter  requirement  is  satisfied  for  the  angles  of 
interest.  The  series  formula  for  the  hypergeometric  function  of  equation  (A-5)  can  be 
deduced  from  the  general  formula  (compare  chapter  15  of  the  reference). 

Fia,b,c,Zo)~2^  -  ,  (^.6) 

n=0 

where  (a)n,  ib)„,  and  (c)„  are  Pochhammer  symbols  defined  by  the  equation 

_r(^  +  n) 


(^)n  = 


(A-7) 


with  r  being  the  gamma  function. 
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Setting 


a  =  -u  , 

(a) 

b  =  u  +  l  , 

(b) 

c=  1  , 

(c) 

and  using  the  properties  of  the  gamma  function  gives 


(a)n  =  (-u)(-u  +  \)(-u  +  2) ...  (-M  +  /i  -  1)  ,  (a) 

ib)n  =  (u+l)iu  +  2) ...  (u  +  n)  .  (b) 

The  resulting  series  expression  for  Pui^)  is 


I’M  =  ‘  £ 

n=l 


ia\(b)„  (I-zq)" 
(«!)'  2" 


(A-8) 


(A-9) 


(A-10) 


When  equation  (A-10)  and  its  counterpart  for  Pu(-2q)  are  inserted  in  equation 
(A-4),  we  have  a  polynomial  of  infinite  order  for  the  determination  of  the  roots  a/.  Since 
this  is  impossible  to  deal  with,  it  is  necessary  to  approximate  Pu(zo)  and  Pu(-zo)  by  a  finite 
number  of  terms. 

In  all  cases  equation  (A-10)  is  a  well-behaved  representation  of  Pui^o)-  Figures 
A-1  and  A-2  show  the  behavior  of  Fu(zo)  and  Pui-^0%  respectively,  for  selected  values  of 
ZQ  =  cos  00  calculated  from  equation  (A-10)  using  128  terms  in  the  summation  and 
including  double  precision. 

The  need  for  double  precision  in  the  computation  of  Pui^o)  and  Pu(-zq)  arises  from 
the  oscillatory  nature  of  the  individual  terms  in  the  series,  which  can  become  extremely 
large  for  correspondingly  large  values  of  u.  Thus,  in  the  absence  of  double  precision,  we 
would  be  faced  with  large  round-off  errors  resulting  from  the  subtraction  of  sequences  of 
two  extremely  large  numbers.  These  uncertainties  would  be  pronounced  in  the 
determination  of  the  roots  of  equation  (A-4).  The  use  of  double  precision  combined  with 
128  terms  in  the  series  of  Pui^o)  and  Fu(-zo)  appears  to  circumvent  the  aforementioned 
problem. 

For  very  large  values  of  u,  even  128  terms  may  not  be  sufficient  to  accurately 
compute  Pu{i-2o).  Fortunately,  in  this  case  the  roots  of  equation  (A-4)  may  be  determined 
using  the  asymptotic  f'^rms  of  Pu{±zq).  The  asymptotic  expressions  [2]  for  the  Legendre 
functions  also  provide  an  excellent  starting  guess  for  the  ZBRAC  and  RTBIS  root  finding 
algorithms  [3]  which  are  used  to  solve  equation  (A-6). 
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-  THETAO  •  16  DEQ  -  THETAO  •  30  DEQ  -  THETAO  -  46  DEQ 

THETAO  -  60  DEQ  -  THETAO  •  76  DEQ 

Figure  A-1.  Pu(2q)  as  a  function  of  u  with  zg  as  a  parameter. 
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-  THETAO  •  16  DEQ  -  THETAO  -  30  DEQ  -  THETAO  ■  46  DEQ 

THETAO  •  80  DEQ  -  THETAO  -  76  DEQ 

Figure  A-2.  Pui-^o)  as  a  function  of  u  with  zg  as  a  parameter. 
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Several  asymptotic  expressions  were  examined,  and  the  one  which  seems  most 
appropriate  for  this  investigation  is  that  rendered  by  Magnus  and  Oberhettinger  [2].  The 
asymptotic  expression  for  the  Legendre  function  of  order  u  is  denoted  by  P„(cos  &)  and  is 
given  by 

P^(cose)=  /  _^  cosHu  +h  d-^;ie=6=n-e,e>0,\u.\»^) .  (A-11) 
mt.  sin  0  L  2  4 j  I  e 

Using  equation  (A- 1 1)  in  equation  (A-4)  and  recalling  that 

-ZQ  =  cos(;r-  0o)  (A- 12) 

gives  the  following  equation  for  the  asymptotic  roots  i7. : 

cos[(m.  +  ^)  00-  -cos[(w:+  ~)in  ■  (A-13) 

Using  the  formula 

cos  a  -  cos  ft  =  -2  sin  ^  {a  +ft)]  sin  [-^  (a  -ft)]  (A-14) 


and  ietting 


d  =  0.00 


(a) 


b  =  vfK-eQ)-^  ,  (5)  (A-15) 


reduces  equation  (A-13)  to  the  form 

sinCw  y)  sin[o.(-^-0o)]=O  .  (A-16) 

Since  we  are  looking  for  solutions  where  u.  is  not  an  integer,  the  solution  of  equation 
(A-16)  is  given  by 

v.(^-eo)=m7i:  ,  (A-17) 

where  m  is  an  integer.  Letting 
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(A-18) 


gives 


-^0 


1 

2 


(A-19) 


As  observed  from  equation  (A-19),  every  angle  do  generates  an  infinite  set  of  roots 
with  magnitudes  increasing  with  the  index  m.  Table  A-1  compares  the  asymptotic  roots 
given  by  equation  (A-19)  and  the  exact  roots  determined  from  the  numerical  solution  of 
equation  (A-4). 

The  asymptotic  values  appear  to  provide  a  good  approximation  at  the  smaller  values 
of  Ui  with  excellent  agreement  occurring  at  the  larger  ones.  This  is  consistent  with  the 
theoretical  expectations. 


Table  A-1.  Comparison  between  exact  and  asymptotic  roots 


do  (degrees) 

Root  no. 

Asymptotic  value,  u. 

Exact  value  Uj 

15.0 

1 

1.90000 

1.80156 

15.0 

2 

4.30000 

4.25625 

15.0 

3 

6.70000 

6.67405 

15.0 

4 

9.10000 

9.04312 

15.0 

5 

11.5000 

11.4272 

30.0 

1 

2.50000 

2.44524 

30.0 

2 

5.50000 

5.45228 

30.0 

3 

8.50000 

8.45644 

30.0 

4 

11.5000 

1 1 .4999 

30.0 

5 

14.5000 

14.4999 

45.0 

1 

3.50000 

3.45397 

45.0 

2 

7.50000 

7.48620 

45.0 

3 

11.5000 

11.4812 

45.0 

4 

15.5000 

15.4990 

45.0 

5 

19.5000 

19.4997 

60.0 

1 

5.50000 

5.48263 

60.0 

2 

11.5000 

11.4999 

60.0 

3 

17.5000 

17.4993 

60.0 

4 

23.5000 

23.5000 

60.0 

5 

29.5000 

29.5001 

75.0 

1 

11.5000 

11.4870 

75.0 

2 

23.5000 

23.4998 

75.0 

3 

35.5000 

35.5005 

75.0 

4 

47.5000 

47.5000 

75.0 

5 

59.5000 

59.4999 
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It  should  be  noted  in  passing  that  Schelkunoff  also  obtained  estimates  for  the 
asymptotic  roots  [4]  which  agree  with  equation  (A- 19)  only  in  the  limit  where  » 
1/2.  He  did  not  attempt  to  generate  numerical  solutions  relevant  to  our  range  of  interest 

After  the  roots  are  determined  it  is  also  necessary  to  calculate  the  angular 
dependence  of  the  function 

(cos  ^  =  4-  [^u.(cos  6)  (-cos  (A-20) 

which  from  equation  (A-4)  is  observed  to  satisfy  the  condition 

(cos  Bq)  =  [cos(tc  -  0o)]  =  (-cos  6q)=0  (A-21) 

i  i  i 

Figures  A-3  to  A-7  show  plots  of  Tu^(z  =  cos  0),  Ptt,(z  =  cos  6),  and 
=  “COS  6)  as  a  function  of  6  for  the  first  five  roots  when  Qq  =  30°.  Using  equation 
(A- 10),  we  can  compute  the  angular  dependence  of  r„j(z)  from  the  expression 
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7’„(cos^=±^^^^^^4j[(l-cos0)'‘-(l+cose)"]  ,  (A-22) 

n=i  (n\)  2 

where  the  (a)n  and  (b)n  are  calculated  from  equation  (A-9)  with  the  m/’s  given  in  table  A-1 
for  6b  =  30°. 

It  is  also  of  interest  to  examine  the  sensitivity  of  the  behavior  of  Tuiicos  6)  to  the 
choice  of  the  asymptotic  versus  exact  root  as  a  function  of  6b.  Figures  A-8  to  A- 12  show 

the  comparison  between  r„-  and  for  the  first  five  modes  when  =  30°.  As  expected, 

the  differences  become  smaller  as  m  increases  from  1  to  5. 

Figure  A- 13  shows  the  comparison  between  the  Tu/s  for  the  smallest  root,  which 
occurs  in  the  6b  =  15°,  m  =  1  case.  For  this  situation  it  is  observed  that  the  difference  can 
become  quite  large,  indicating  that  the  asymptotic  approximation  is  not  especially  good.  On 
the  other  hand,  figure  A- 14  shows  the  Oq  =  75°,  m  =  1  case,  which  corresponds  to  a 
relatively  large  root  of  11.50  (compare  table  A-1).  As  observed  from  figure  A-14  the 
differences  in  this  case  are  imperceptible. 
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0  20  40  60  80  100 

0 

- P(2)  - P(-Z)  - T(Z)-.5  (P(Z)-P(-Z)) 


Figure  A-3.  Tu(cob  0),  Pu(cos  0),  and  Pu(-cos  0)  as  a  function  of  0  for  0qs  30° 

and  m  s  1. 


- P(Z)  —  P(-Z)  - T(ZK6  (P(Z)-P(-Z)) 

Figure  A-4.  Tuicos  0),  P|/(cos  0),  and  P{/(-cos  0)  as  a  function  of  0  for  0q  s  30° 

and  m  =  2. 
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0  20 


40  60  80  100 

9 


P(Z)  -  -  -  P(-Z)  - T(Z)-.5  (P(Z)-P(-Z)) 


Figure  A-5.  Tuicos  9),  Pu(cos  9),  and  Po{-cos  9)  as  a  function  of  9  for  0o=  30“ 

and  m  s  3. 


- P(Z)  -  -  P(-Z)  - T(Z)-.6  (P(Z)-P(-Z)) 


Figure  A-6.  Tuicos  6),  Puicos  9),  and  Pu{-cos  9)  as  a  function  of  9  for  dgs  30° 

and  m  s  4. 
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T^cos  6)  3  Tuicos  0),  Puicos  eh  Pul-cos  B) 


0  20  40  60  80  100 

e 


- P(Z)  -  -  -  P(-Z)  - T(Z)-.5  (P(Z)-P(-Z)) 


gure  A-7. 


Tuicos  0),  Pulcos  6),  and  Pul-cos  6)  as  a  function  of  6  for  Oq  s  30° 
and  m  a  5. 


0  20  40  60  80  100 

9 

- ASYMPTOTIC  -  -  EXACT 

Figure  A-8.  Comparison  between  asymptotic  and  exact  Tu  lor  ©gr  30° 

and  m  s  1. 
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TiXcose)  Tj/COSG) 


- ASYMPTOTIC  -  -  -  EXACT 

Figure  A>9.  Comparison  between  asymptotic  and  exact  Tu  for  Oq  =  30° 

and  m  s  2. 


- ASYMPTOTIC  -  EXACT 

Figure  A>10.  Comparison  between  asymptotic  and  exact  Tu  for  0q  s  30° 

and  m  s  3. 
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.3 


ASYMPTOTIC  -  -  EXACT 


Figure  A-11.  Comparison  between  asymptotic  and  exact  7^  for  do=  30° 

and  m  s  4. 


0  20  40  60  80  100 


9 

- ASYMPTOTIC  -  EXACT 

Figure  A-12.  Comparison  between  asymptotic  and  exact  Tu  for  dgs  30° 

and  m  s  5. 
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7*1/008  0)  Tu(COS  0) 


0.2 


0.16 

0.1 

0.06 

0 

■0.06 

0  20  40  60  80  100 

6 

- ASYMPTOTIC  -  -  EXACT 

Figure  A-13.  Comparison  between  asymptotic  and  exact  for  9q  =  15° 

and  m  =  1. 


0  20  40  60  80  100 

e 

- ASYMPTOTIC  ---EXACT 

Figure  A-14.  Comparison  between  asymptotic  and  exact  Ty  for  0o  s  75° 

and  m  s  1. 
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